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The friction coefficient of a particle can depend on its position as it does when the particle is near a 
wall. We formulate the dynamics of particles with such state-dependent friction coefficients in terms 
of a general Langevin equation with multiplicative noise, whose evaluation requires the introduction 
of specific rules. Two common conventions, the Ito and the Stratonovich, provide alternative rules 
for evaluation of the noise, but other conventions are possible. We show the requirement that a 
particle's distribution function approach the Boltzmann distribution at long times dictates that 
a drift term must be added to the Langevin equation. This drift term is proportional to the 
derivative of the diffusion coefficient times a factor that depends on the convention used to define 
the multiplicative noise. We explore the consequences of this result in a number examples with 
spatially varying diffusion coefficients. We also derive path integral representations for arbitrary 
interpretation of the noise, and use it in a perturbative study of correlations in a simple system. 

PACS numbers: 05.40.-a 



I. INTRODUCTION 

Brownian motion provides a paradigm for exploring 
the dynamics of nonequilibrium systems, especially those 
that are not driven too far from equilibrium [J [2j, H, H| • 
In particular, the Langevin formulation of Brownian mo- 
tion finds applications that go beyond its original purpose 
of describing a micron-sized particle diffusing in water. 
It has been extended to treat problems in dynamics of 
critical phenomena in glassy systems @, and even 
in evolutionary biology 0]. Brownian motion is impor- 
tant for soft-matter and biological systems because they 
are particularly prone to thermal fluctuations [l|, Q , and 
Langevin theory is an important tool for describing their 
properties, such as the dynamics of molecular motors Q 
and the viscoelasticity of a polymer network [lfj . 

In most applications, the diffusion coefficient is as- 
sumed to be independent of the state of the system. Yet, 
there are many soft-matter systems in which the diffusion 
coefficient is state dependent. A simple example of such a 
system is a particle in suspension near a wall: its friction 
coefficient, and hence its diffusion coefficient, depends be- 
cause of hydrodynamic interactions on its distance from 



the wall a phenomenon that affects interpretation 



of certain single-molecule force-extension measurements 
[l2| and that plays a crucial role in experimental verifica- 
tion of the fluctuation theorem in a dilute colloidal sus- 
pension near a wall [l3j]. Similarly, the mutual diffusion 
coefficient of two particles in suspension depends on their 
separation p3 |. Other examples with state-dependent 
diffusion include a particle diffusing in a reversible chemi- 
cal polymer gel [15| and the dynamics of fluid membranes 
(l6| . In spite of the recent advances in digital imaging 
methods to probe equilibrium properties of soft matter 
[l7| . there have been relatively few experimental studies 



of the dynamical properties of a physical system in which 
the diffusion coefficient is state dependent. This is clearly 
an area for further experimental exploration. Although 
the mathematical problem of how to treat systems with 
state-dependent diffusion has been studied for some time 
0, [l(| H, O H3| ; the results of these studies have not 
been collected in one place to provide a clear and concise 
guide to both theorists and experimentalists who might 
use them. 

In this largely expository paper, we develop a Langevin 
theory and its associated path integral representation for 
systems with state-dependent diffusion and explore its 
use in systems ofphysical interest. In accord with previ- 
ous treatments [J 0, E3, EH, El, QjJ , we show that a 
position-dependent diffusion coefficient leads naturally to 
multiplicative noise. This noise is the product of a state- 
dependent prefactor proportional to the square root of 
the diffusion coefficient and a state-independent depen- 
dent Gaussian white noise function, and it is meaningless 
without a prescription for the temporal order in which 
the two terms are evaluated. There are two common 
prescriptions or conventions for dealing with multiplica- 
tive noise: the Ito convention in which the prefactor is 
evaluated before the Gaussian noise and the Stratonovich 
convention which results when the delta-correlated white 
noise is obtained as a limit of a noise with a nonzero 
correlation time 0, S 0|- There are, however, other 
conventions as we will discuss. Using general thermo- 
dynamic arguments, we show that in order for Boltz- 
mann equilibrium to be reached a drift term propor- 
tional to the derivative of the diffusion coefficient times 
a factor depending on the convention for the evaluation 
of multiplicative noise must be added to the Langevin 
equation. Though this drift term has been noted before 
we have found only one (recent) reference 
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1 1 01 ] that specifically associate the form of the drift term 
with the convention for evaluating multiplicative noise. 
On the other hand, others claim that it is the choice of 
the convention that is dictated by physics [ll[. In par- 
ticular, the authors of Ref. without allowing for the 
possibility of the drift term, argued that neither Ito nor 
Stratonovich convention properly describes the dynamics 
of a Brownian particle with a spatially varying friction 
coefficient, but a third convention - what the authors 
called the isothermal convention, does. Incidentally, for 
this third convention, the drift term in our formalism 
vanishes. Therefore, the necessity of the drift term for 
enforcing thermal equilibrium is not widely known, and 
it is often incorrectly ignored [T3 |. Here, we aim to pro- 
vide a clear exposition for clarifying the technical issues 
that might have been a source of confusion in the litera- 
ture. 

This paper is organized as follows: in Sec. [TTJ we first 
review the case of a uniform diffusion coefficient and ex- 
tend it to the case of spatially varying diffusion coef- 
ficient. We discuss in depth the stochastic interpreta- 
tion issues associated with multiplicative noise, we de- 
rive the Fokker-Planck equation, and we show that de- 
pending on the stochastic interpretation, an additional 
drift term must be added to the standard friction term 
in order for the system to relax to equilibrium. We also 
discuss how measurements of the eigenvalues and eigen- 
functions of the probability that a particle is at position 
x' at time t + St given that it was at position x at time 
f can be used to obtain information about whether the 
diffusion coefficient is state-dependent or not. In Sec. lIIIl 
we present some exactly solvable toy models that clearly 
illustrate the consequences of spatially varying diffusion 
and suggest some experimental techniques which may 
elucidate its role in colloidal tracking experiments. We 
also give numerical confirmation that the extra drift term 
is needed to produce equilibrium distribution. In Sec. 
TVl we derive and discuss the path integral formulation 
for a Langevin equation with a multiplicative noise, cor- 
relation functions, and perturbation theory. In Sec. [Vj 
we briefly summarize the results for multicomponent sys- 
tems. Technical details are presented in the Appendices. 



II. FORMALISM IN 1-D 

A. A review for the case of a uniform diffusion 
coefficient 

Let us first briefly review the simplest case in which a 
Brownian particle diffuses in space with a uniform diffu- 
sion constant [2j. In the Langevin formulation of Brow- 
nian motion, the stochastic equation of motion for the 
particle's position Q is 



where x denotes the position, T is the dissipative co- 
efficient (inverse mobility), Ti. is the Hamiltonian, and 
g 77(f) models the stochastic force arising from the rapid 
collisions of the water molecules with the particle. The 
strength of this force is set by g, and 77(f) is a Gaussian 
white noise with zero mean, (77(f)) = and variance, 
(r)(t)r)(t')) = 6(t — f'), delta-correlated in time. The first 
term on the right hand side of Eq. (|2.ip describes a dissi- 
pative process. Thus, Eq. (|2.ip can be viewed as a balanc- 
ing equation in which the first term drains the energy of 
the particle while the random noise pumps it back. Equa- 
tion (|2.1[) neglects an inertial term that is only important 
at short times, typically less than 10" 7 s in soft-matter 
systems [§]. Thus, Eq. ()2.1|) tacitly assumes that there is 
a separation of time scales in which the time scale of the 
fast processes reflecting microscopic degrees of freedom is 
much shorter than the typical time scale for the random 
variable x(t). Hence, the white noise assumption in Eq. 
(EU). 

The Fokker-Planck equation 0,0, 0], 



d t P(x,t)=d x 



&H 

dx 



1 



9 2 d x 



P(x,t), 



(2.2) 



d t x 



an 

dx 



(2.1) 



which can be derived for the Langevin equation, for the 
probability density P(x, f) that a particle is at position 
x at time f provides an alternative to Langevin equation 
for describing the motion of Brownian particles. It is 
easy to see that Eq. (|2.2[) has a steady state solution 
P s (x) oc exp[— 7i]. If a particle is in equilibrium with 
a heat bath at temperature T, then P 8 {x) oc exp[— /37i] 
from which we conclude that g 2 — 2TkBT. If H = 0, 
Eq. (|2.2|) reduces to a diffusion equation with diffusion 
constant D = g 2 /2. Hence, for systems in equilibrium at 
temperature T, the diffusion constant obeys the Einstein 
relation D = k B TT. 



Extension to the case of state-dependent 
diffusion coefficient 



How must the Langevin equation for a Brownian par- 
ticle be modified when the friction coefficient T de- 
pends on position xit), i.e., when T depends on the 
state of the system.? Though it is generally understood 
[lUilUlHilMlH that an x-dependent T leads to an in- 
dependent g and thus to multiplicative noise g [x(t)} 77(f), 
it is less well known that the requirements of long-time 
thermal equilibrium require an additional specific mod- 
ification to the Langevin equation - the addition of a 
convention-dependent drift term. Though there are dis- 
cussions in the literature of this drift term [Tol. [ill [lH [20| . 
they are not very detailed, and they generally treat only a 
specific convention for dealing with multiplicative noise. 
Here we show that constraints of equilibrium require a 
unique drift term with each noise convention and resolve 
any ambiguities [22| arising from the fact that multiplica- 
tive noise can be interpreted in many ways. 
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Using the argument that the stochastic force is bal- 
anced by the dissipative term as in the case of a uniform 
dissipative coefficient above, we may reasonably postu- 
late a Langevin equation, which trivially generalizes Eq. 
(|2.f p to the case of spatially varying dissipative coeffi- 
cient, to take the following form: 

d t x = -T(x)^+g(x)r 1 (t), (2.3) 

where g(x) = ^2k B TT(x). But we must first confront 
the issue of interpreting the multiplicative noise g(x)rj(t), 
which by itself is not defined [2j, |2l|. This is because 
the stochastic nature of rj(t) which in general consists 
of a series of delta-function spikes of random sign. The 
value of g[x(t)]rj(t) depends on whether g[x(i)] is to be 
evaluated before a given spike, after it, or according to 
some other rule. It turns out, as we will show shortly, 
that this naive generalization of Eq. (|2 . 1 [) to Eq. (|2.3[) is 
only valid for a particular interpretation of the noise. 

There are a number of approaches to assigning mean- 
ing to the multiplicative noise, but they all boil down to 
providing rules for the evaluation of the integral 

/t+At 
dsg[x(s)]r,(s), (2.4) 

in the limit of small At. If g(x) and 77(5) are both contin- 
uous functions, this integral could, for arbitrary At, be 
expressed via the first integral mean-value theorem as 

/t+At 
ds r){s) (2.5) 

where U is a uniquely determined time in the interval 
[t,t + At]. In the limit of small At, this expression, Eq. 
(|2.5p , reduces trivially to g[x(t)]rj(t) At, to lowest order in 
At. The noise r](s) is, however, not continuous and Eq. 
(12. 5p with a uniquely determined time does not apply. 
One can, however, use Eq. (|2.5p to motivate a definition 
of J(t, At) for a stochastic rj(s). There are two commonly 
used conventions for defining 3{t, At): the Stratonovich 
convention 

/t+At 
dsri(s), (2.6) 

in which g[x(t)] is evaluated at the midpoint of the inter- 
val [x(t),x(t + At)] and the Ito convention 

rt+At 

Ji(t) = g[x(t)] J dsr](s), (2.7) 

in which g[x(t)] is evaluated before any noise in the inter- 
val (t, t+At) occurs. We will use a generalized definition: 

/t+At 
dsri(s), 

(2.8) 



which is parameterized by a continuous variable a € 
[0,1], that reduces to the Ito convention when a = 0, 
to the Stratonovich convention when a = 1/2, and to 
the isothermal convention of Ref. [llj when a = 1. 

We note in passing that in the mathematics commu- 
nity, the Ito calculus is most commonly used. Perhaps, 
this is because of the conceptual simplicity arising from 
the property that the noise increment J^ +At dsij(s) and 
x(t) are statistically independent as implied in Eq. (|2.7[) . 
i.e. (g(x)r](t)) = [231] . On the other hand, in the 
physics community, the Stratonovich interpretation is fa- 
vored. In addition to the advantage that it gives rise to 
the ordinary rules of calculus, the Stratonovich conven- 
tion also has a deeper physical origin. Since the noise 
term in Eq. (|2.3p models, in a coarse-grained sense, the 
effects of microscopic degrees of freedom that have finite 
(albeit short) correlation times, this term should be phys- 
ically interpreted as the limit in which these correlation 
times go to zero. By the Wong-Zakai theorem, this limit 
corresponds to a white noise that must be interpreted us- 
ing the Stratonovich convention [23j | . However, Eq. (|2.3p 
does not provide a correct description for systems in con- 
stact with a thermal bath at temperature T for cither 
interpretation: their associated Fokker-Planck equations 
do not have long-time thermal-equilibrium solutions. 

To return to our main discussion, it is clear that 
J a (t,At) depends on the value of a. Integration of Eq. 
(E3D yields 

Ax{t + At) = x(t + At) - x(t) = J a (t, At), (2.9) 

when Tt = 0. The integral J* +At dsr](s) is statistically 
of the order of \/Ai, implying Ax(t + At) is also of the 
order of y/Al. Thus, ax(t + At) + (1 - a)x(t) = x(t) + 
aAx(t + At) has a term of order y/At proportional to 
a, and the order At term in jf a (t, At) depends on a. 
An alternative approach to defining J(t, At) is simply 
to expand g[x(s)] in the integrand as g[a;(s)] = <7[cc(t)] + 
[x(s) — x(t)] g'[x(t)] + ■ ■ ■ . In this approach, which we 
outline in Appendix [Al ambiguities in the interpretation 
of J(t, At) are resolved by specifying the value of the 
Heaviside unit step function, 6(t) at t = 0. Setting 9(0) — 
a is equivalent to using Eq. (|2.8p for J(t, At). 

The stochastic integral J a (t, At) depends on our con- 
vention for evaluating it, i.e. on a. Thus, different val- 
ues of a define different dynamics. But the requirements 
of thermal equilibrium should imply a unique dynamics. 
What is missing? To resolve this dilemma, we consider 
the general stochastic equation 

d t x = f(x)+g(x)r 1 (t), (2.10) 

where 

f( x ) = -T(x)^ + f 1 (x), (2.11) 

in which we leave f\(x) unspecified for the moment. Eq. 
(I2.10p is easily integrated using the rules we just outlined 
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to yield 
x(t + At) - x{t) 



t+At 



da{f[x(a)]+g[x(s)] V (s)} 



f[x(t) +aAx]At 



+ g[x{t) + aAx] 



t+At 



dsr)(s)l2.12) 



from which we obtain, to the first order in At, 

(Ax) = f(xo)At + ag(x )g'(x )At, (2.13) 
((Ax) 2 ) = .g 2 (x )At, (2.14) 

where we set x(t) = Xq. Thus, there is a stochastic 
contribution, agg' At, to (Ax) arising from the x depen- 
dence of g and depending on the convention for evalu- 
ating J(t, At). In equilibrium, (Ax) should be indepen- 
dent of a. Thus, it is apparently necessary to include a 
contribution to /(x) depending on a. 



C. Derivation of the Fokker-Planck Equation and 
Equilibrium Conditions 

To determine the appropriate form of /(x) and g(x) 
to describe equilibrium systems with a spatially vary- 



ing friction coefficient T(x), we derive the Fokker-Planck 
equation for the probability density P(x, t). The Fokker- 
Planck equation is most easily derived using the identity 



P(x, t + At) = J dx P(x, t + At\x t)P(x t), (2.15) 



where P(x, t + At\xo t) is the conditional probability dis- 
tribution of x at time t + At given that it was xo at time 
t. It is defined by 



P(x, t + At\x a t) = (S[x - x(t + At)}} 



x ,t 



(2.16) 



where the average is over the random noise g(s) and x(t+ 
At) is determined by Eq. (|2.12|) with x(t) = x . Taylor 
expanding the conditional probability around xq yields 



P(x,t + At\x t) 



d 

S(x - X ) - {Ax} — S(x - Xq) 

ox 

^((Ax) 2 )^^-x ) + .... 



J 



Then using this in Eq. (|2.15p . we obtain 







d t P(x,t) = ^[-f(x)-ag(x)g'(x)}P(x,t) 



i d 2 

2dx 2 



g 2 {x)P{x,t)} 



0_ 

dx 



P{x,t). 



(2.17) 
(2.18) 



For an equilibrium system, this equation must have a 
steady state solution with the canonical form 



P{x,t) 



-n/{k B T) 



(2.19) 



that is always approached at long times. Such a solution 
is guaranteed if 



ff 2 (z) 
fi( x ) 



2k B TT(x), 

(1 - a)g(x)g'(x) 



2(1 



(2.20) 

a)k B TT'(x).{2.2\) 



Thus, an additional drift term, /i(x), which depends on 
the convention for evaluating Jit, At), must be added to 
the standard friction term, — T(x)d x Tl, in the equation 
for dtx in order for the system to evolve to the Boltz- 
mann distribution at long times, i.e., be consistent with 
thermodynamics. Note that /i(x) is proportional to the 
temperature T, indicating that its origin arises from ran- 
dom fluctuations rather from forces identified with a po- 
tential. It is clear now from Eq. (|2.21[) that if we insist 
on using the Langevin equation in the form of Eq. 
we are forced to take a = 1 [ill. 



It is customary to express the Fokker-Planck equa- 
tion in terms of the diffusion constant rather then the 
friction coefficient. From Eq. (|2.14[) for ((Ax) 2 ), we 
can identify g 2 (x) with the short-time diffusion constant 
D(x) = 2k B TT(x). With this definition of D(x) and 
/i(x) given by Eq. (|2.21[) . the Fokker-Planck equation 
becomes 



d t P(x,t) = ^D(x) 



(3— + — 

dx dx 



P(x,t), (2.22) 



where = l/(fcsT). As required, this equation is inde- 
pendent of a: different conventions now give the same 
equilibrium condition as they should. For a free particle 
diffusing in spatially varying D(x), Tt = and Eq. (|2.22|) 
becomes 

d t P(x, t) = — D(x) — P(x, t). (2.23) 

This implies that the correct generalization of Fick's Law 
for equilibrium systems with a spatially- varying diffusion 
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coefficient is given by 



III. ILLUSTRATIVE EXAMPLES 



J(x,t) = -D(x)d x P(x,t). 



(2.24) 



Historically, the generalization of Fick's law has long been 
debated [24]]. It is commonly acknowledged that Eq. 
(|2.24[) is right even though many derivations to the right 
of side of Eq. (12.24[) seem not to be as transparent as the 
one given above. 



In this section, we consider some exactly solvable toy 
models to illustrate some central ideas presented in the 
last section. In particular, we address the effects of spa- 
tial dependence in the diffusion coefficient and use numer- 
ical solution of the Langevin equation to show that equi- 
librium distribution is obtained only if fi(x), Eq. (|2.21[) . 
is added to the standard friction term. 



D. Experimental probes of D(x) 



One interesting property of Eq. (|2.22|) is that it nec- 
essarily has an eigenstate with eigenvalue zero and 
eigenfunction given by the equilibrium distribution the 
P e „(x) oc e~P HtyX \ This fact is exploited by Crocker et 
al. 17[ to measure directly the interaction between an iso- 
lated pair of colloidal particles. In these experiments, the 
data from tracking the motion of the particles are used 
to compute the conditional probability P(x,t + St\x't), 
which may be viewed as the Green's function to or the 
inverse of the Fokker-Planck equation, Eq. (|2.22[) . The 
equilibrium distribution is then the solution to 



P eq {x) = / dr' P(x, t + St\ x', t)P eq (x'), 



(2.25) 



from which the interaction potential can be constructed 
via U(x) = -k B T dog P eq (x). 

Since the conditional probability contains all the dy- 
namical information of the system, one could in princi- 
ple characterize how the system relaxes to equilibrium 
by extracting the nonzero eigenvalues of Eq. (|2.22| . In 
particular, the Fokker-Planck equation describing a sys- 
tem with a state-dependent diffusion coefficient would 
have eigenvalues and eigenfunctions that are, in general, 
different from those of a system with a uniform diffu- 
sion coefficient, even though the two systems have the 
same Hamiltonian. Thus, in principle, one could extract 
the Hamiltonian from image analysis following the proce- 
dures of Crocker [13], and solve Eq. (|2.22p with uniform 
diffusion constant to obtain a set of eigenvalues (prob- 
ably numerically) and compare it with experimentally 
measured eigenvalues, which can be extracted from the 
measured conditional probability. If they are different, 
then the diffusion coefficient is state dependent, and one 
needs to model the diffusion coefficient to understand the 
dynamical behaviors of the system. We suggest this pro- 
cedure as a possible general method for experimentalists 
to explore the dynamics and measure the state dependent 
friction coefficient in, for example, hydrodynamic inter- 
actions between two spheres [25[ , diffusion of particles in 
polymer solution [26| , and rods in a nematic environment 



A. Diffusion of a particle near a wall 

How does a diffusion coefficient acquire a spatial de- 
pendence? The simplest example is a Brownian particle 
diffusing near a wall located at z = 0. Brenner [28| has 
shown that for z > the diffusion coefficient acquires a 
spatial dependence in which it is zero at the wall, rises 
linearly in z, and approaches a uniform bulk value of Do 
at large z as 



D(z) 



9a 

D ~ 8z 



(3.1) 



Note the long-range component of D(z) in Eq. (|3.ip . 
which reflects the long-ranged nature of the hydrody- 
namic interaction. Recently, it has been pointed out that 
in single molecule experiments, it is crucial to take the 
spatial dependence in the diffusion coefficient properly 
into account [l2T ] . 

Rather than to treat the system with the above D(z), 
we consider a toy model, which correctly describes diffu- 
sion close to a wall, in which D(z) — Az. This diffusion 
coefficient has another experimental realization: diffusion 
of a colloidal particle bounded by two parallel walls, with 
one of the walls slightly tilted ll| . Then, the diffusion 



coefficient acquires a spatial dependence, approximately 
given by D(z) ~ z, for the motion of the particle parallel 
to the walls. In this case, the Fokker-Planck equation 
becomes 



d t P(z,t) = A-^-z-^P(z,t), 

az az 



(3.2) 



which can be solved exactly. Let P(z,t) = 
J2n c ne~ Xnt ip n {z), where A„ are a set of the eigenval- 
ues. With the transformation y = y/z, Eq. (|3.2p can be 
written as 

„ i , a\ 

^„(v) + -^(») + -^n(y) = o, (3.3) 

whose solution is the Bessel function: i/j(y) — Jo(ky), and 
whose eigenvalues form a continuous spectrum given by 
A = Afc 2 /4. The probability distribution as a function of 
time can be written as 



P(z,t)= dkc(k)e~ Ak2t/4 J Q (k^z~). 



(3.4) 
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FIG. 1: Stationary distribution for a particle diffusing in a 
diffusion coefficient D(z) = Az subject to a constant force 
F = —Fz. What is shown here is the numerical simulation 
of the Langevin equation Eq. I|3.10[) : it is of the form e~^ Fz , 
as expected. 



The probability distribution for a particle at z = zq at 
time t = evolves as 

P(z, t) = j- t e -[(^«)/(A01 /o [ 2 ^E/(At)] . (3.5) 

Unlike its counterpart for a uniform diffusion coeffi- 
cient, this probability distribution is non-Gaussian. It 
is straightforward to calculate the moments: 

(z(t)) = z + At 
([z(t)-z ] 2 ) = 2At(z + At). 

These behaviors are very different from those of a con- 
stant diffusion. In particular, the mean-squared dis- 
placement exhibits ballistic behavior. It is interesting 
to observe that the second moment can be written as 
((z(t) - z Q ) 2 ) = 2A(z(t))t. This suggests that m or- 
der to extract the diffusion coefficient for this simple 
problem, we need to know not only the second mo- 
ment ((z(t) — Zo) 2 ), but also the first moment (z(t)). 
Only for short times does the second moment reduce to 
{(z(t) — zq) 2 ) ~ 2 Azot — 2D(zo)t, which is the formula 
commonly used to extract the diffusion coefficient. It 
is clearly incorrect to use this formula for times greater 
than Zo/A. The method we suggested at the end of the 
last section compliment this approach. Note also that 
the (z(t)) ~ t behavior has been measured in Ref. [Til ]. 

If the particle is subject to constant force F, like grav- 
ity, in the — z direction, then the Fokker-Planck equation 
is 



d t P(z,t) = A-^-z 
az 



az 



P(z,t). 



(3.6) 



This problem can also be solved exactly. Let P(z,t) = 
e -PFz ^ c n e~ Xnt ip n (z), we find that the eigenfunctions 
satisfy the Laguerre equation 



xip + (1 — x)ip + 



X, 



Af3F 



^ = 0, 



(3.7) 



with eigenvalues A„ = nAfiF. The eigenvalue spectrum 
is discrete rather than continuous as it is in the case of 
a constant diffusion coefficient. If the particle is initially 
at zq, the distribution evolves as 



P(z,t) 



(3F 



I _ e ~A0Ft 
X I 



exp - 



z e 



-A/3F* 



1 _ g-A/3-Ft 



20Fy/ 



zoze 



-A0Ft 



1 



-A0Ft 



(3.8) 



Note that at t — > oo, this distribution reaches the equilib- 
rium distribution P eq ~ e~@ Fz . The first two moments 
of z(t) are 

(*(*)> 



z e 



-AfiFt 



1 



1 



-A0Ft\ 



(z 2 (t)) = z 2 e- 2 ^ Ft 
2 



+ 



(f3F)- 



-(1 



(3F 

4z e- A ^ Ft 
I 

f3F 

e ~A0Ft^2^ 



(1-e 



-A0Ft} 



(3.9) 



Note that at long time (z(t)) = kgT/F as thermal equi- 
librium dictates. 

We numerically solve the Langevin equation corre- 
sponding this problem 



d t z{t) = A + A(3Fz + V2Azr)(t), 



(3.10) 



in the Ito convention [29(. Note that the first term in 
the right-hand side arises from the additional drift, f\{x) 
given by Eq. (|2.2ip . The result for the stationary distri- 
bution is plotted in Fig. [T] Obviously, it agrees with the 
equilibrium distribution P eq ~ e~@ Fz . 



B. Diffusion of a particle bounded by two parallel 
walls 



Next, we consider the diffusion of a particle bounded 
by two walls, which was studied experimentally in Ref. 
[30l ] and more recently in Ref. [3l|. We approximate the 
spatially varying diffusion coefficient of this system by 
D(x) = D [1- {x/L) 2 }. The resulting Fokker-Planck 
equation is 

d t P(x, t) = D JL[l- (x/L) 2 } |- P(x, t), (3.11) 

with boundary conditions that particles cannot penetrate 
the walls, i.e., that the flux at both walls be zero: 

J(x, t) = D \l- (x/L) 2 ] ^- P(x, t) = at x = ±L. 
L J ox 

(3.12) 

Again the solution to this problem differs considerably 
from that with a spatially uniform diffusion coefficient. 
The spectrum is discrete rather than continuous with 
eigenvalues A„ = n(n+ 1)D/L 2 (n = 1, 2, ...) and the as- 
sociated eigenfunctions are Legendre polynomials rather 
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FIG. 2: Stationary distribution for a particle diffusing 
between two walls with a diffusion coefficient D(x) — 
Do [l — (t/L) 2 ]. The solid line is the numerical simulation 
of the Langevin equation in Eq. (|3.14p and the dotted line 
is the numerical simulation of the Langevin equation without 
the extra gg' term. 



than linear combinations of plane waves. The first two 
moments of x(t) are 



(x(t)) = x e- 2D ^ L2 



(x(t) 



2 \ = xle~ 6D " t/L 



L 2 



l- e - 6Dot/L2 )(3.13) 



These moments again are different from the case in which 
the diffusion is uniform. 

We performed numerical simulation of the Langevin 
equation 

d t x = ~2D Q x/L 2 + V2A)(1 - x 2 /L 2 ) r](t), (3.14) 

where the first term arises from the gg' term. In Fig. [2l 
we plot the long time distribution (solid line) which is 
uniform as it should be. We also show the numerical re- 
sults for the case in which we did not add the gg' (dotted 
line). Clearly, we get the wrong answer if we do not add 
the gg' term. 



C. Diffusion constant: D(x) = Do(l +7 a; 2 ) 

As a final example, let us consider a free particle dif- 
fusing with D(x) = Do (1+7 a; 2 ) in the bulk. The Fokker- 
Planck equation is given by 

d t P(x,t) = D — [1 + 7 a; 2 ] g^P(x,t). (3.15) 

Multiplying both sides by x 2 and integrating, we find 

d t (x 2 (t)) =2D + 6A)7 (x 2 (t)), (3.16) 
whose solution is 

(x 2 (t)) = -!-[exp(6A,7*)-l]- ( 3 - 17 ) 

37 



Thus, the second moment grows exponentially with time; 
this peculiar behavior illustrates the dramatic effects of 
the noise in problems with a spatial dependent diffusion 
coefficient. 



IV. PATH-INTEGRAL FORMULATION 

Path-integral formalisms provide an alternative to the 
Fokker-Planck and Langevin equations for the descrip- 
tion of stochastic dynamics. They have the advantage 
that well-established perturbative and non-perturbative 
field-theoretic techniques [H, [HJ can be used to calcu- 
late the effects of nonlinearities. They also provide a 
convenient treatment of correlation and response func- 
tions. The path integral for a state-dependent dissi- 
pation coefficient has been derived previously either in 
the Stratonovich or Ito convention 0, [2^, [H, [Hj]. In 
this section, we derive the path-integral for the general 
a convention, use it along with the detailed balance, 
a condition that any thermal systems must satisfy, to 
shed further insight into the additional drift term de- 
rived in Sec. Ill CI We also discuss equilibrium correla- 
tion and response functions and prove the Fluctuation- 
dissipative theorem for state-dependent diffusion coeffi- 
cient. We then set up perturbation theory for systems 
with a coordinate-dependent friction coefficient. 

The path-integral is based on the statistics of a path 
x(t). We discretize the path into segment Xi = xfa) 
with to < t\ < ■ ■ ■ < t TV— 1 < tf and At = 
t n — t n —i small. The joint probability distribution, 
P(xjv*jV) %N-i *JV-ii • • • i x x to) that x{t) takes on 
values of x\ at time t\, xi at time t^ and so on, given 
that it has value of Xq at time to, is then 

P{xf t f ;x N -i tff-i', ...;xi h\x 1 ) 

= (S[x N - cf)(t N ;xo,to)} . . .6[xi - <f>(tt',x ,t )] ), 

where the average is taken with respect to the noise and 
(j>(ti] Xi-i, ti—i) is the solution to the Langevin equation, 
Eq. (|2.10[) . for x{ti) given that x(U-i) = Xi-\. Since the 
noise in Eq. (|2. L0[) is delta-correlated in time, the noise in 
different time intervals is not correlated, and X4 depends 
only on a;,_i. Thus, we can write 

N 

P(x N t N ; ...;xi tt\x to) = [l(S[xi - <f>(ti\Xi-i,ti-i)]). 

i=i 

The function 

P{x i ti\x i -iti- 1 ) = (S[xi - (j>(ti\Xi- 1,^-1 )]) (4.1) 

gives the conditional probability that the random vari- 
able x{t) has the value Xi at time f, given that it had 
a value Xi-\ at Using Eq. (14. ip and the identity 
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P(xt) = J dx' P(xt;x' t'), it is easy to see that 



P(Xi U\xi- 2 k-2) = J dXi^i P(Xi ti\Xi-l ti-i) 

x P{xi-i U-i\xi-2 U- 2 ), (4.2) 
P(xi U) = / dx t -i P{xi U\xi-i t i - 1 )P(x i - 1 ^-0.3) 



Equation (14. 2p is the Chapman-Kolmogorov equation, 
which defines a Markov process [2|, while Eq. (|4.3[) is 
just an identity, true for all stochastic processes. Note 
that a Markov process is completely specified if we know 
P(xi U) and P(xi ti\xi-\ U-i), but they are not arbitrary 
because they are linked through Eqs. (|4.2[) and (|4.3[) . Us- 
ing Eq. (|4.2[) . the conditional probability for the particle 
to go from xq at time to to Xf at time i/ is 



P(xj tf\x a t ) = J dxN-i...J dxi x (4.4) 

(S[x f - ^(t/;a;jv-i,iiv-i)]) ■ • • (<%i - 0(ti;xo>*o)]}- 

This is the basic construct for the path integral. First, we 
have to evaluate (S[xi — 4>{ti\ x^i^U-i)]) . We discretize 
Eq. (j2~T0|) as follows: 

x, = x^ + At f, + g t f dt'r](t'), (4.5) 

where At = ti-ti-x, U = f[axi + (l-a) Xj-i], and = 
g[axi + (1 — a)xi_i]. We next introduce the function 
h(xi,Xi-i): 



h{xi,Xi-i) = 



dt'r)(t'), (4.6) 



which vanishes when xt is the unique solution to the Eq. 
(|4~5]) . (f>(ti;Xi-i,U-i), i.e., h[ct>(ti;xi-i,ti-i),Xi-i] = 0. 
Using the property of the delta function 



5[h(xi,Xi-i)] = 



dh 



dxi 



=<t>(U 



and noting that 



depends only on Xi and 



Xi-i, which are set by the delta function and not ex- 
plicitly on the noise, we have 



(S[h(xi,Xi_i)]) 



dh 



dxi 



S[xi - cj){ti)} 



dh 



dxi 



Xi = <j>(ti) 



(6[ Xi -</>&)]), 



since for any function q[(j>(ti)}, (q[4>(ti)]S[xi — (j>(ti)]) — 
q{xi)(8[xi — (f>(ti)]). We can, therefore, write the condi- 
tional probability as 

P(x i ti\x i - 1 t i - 1 ) = (S[xi - (p(U; Xi-i, t»-i)]) 
dh{xi,Xi-i) 

{d[h(xi,Xi-i)\), 



with 

dh _ 1 

dxi gi 



dxi 



1-aAtfi-a^ ( Xi - Xl - X -Atfi) 

9* 



where ' denotes the derivative. The average over noise 
can be easily done with the aid of the Fourier represen- 
tation of the delta function: 

<W(*i,*«)]> = f^e +1 |(— At « 

J 27T 



x < e 



2tt 6 ' 



where we have made use of the fact that i dt' r)(t') 
is a zero-mean Gaussian random variable with variance 
At. Putting these results together, we can express 
P{xiti\Xi-xti-i) as 



P^Xiti \xi—iti— 1 1 



dki 
2ngi 



1-aAtfl-a^- (xi - - At f l 
9i 



(4.7) 



Next, in order to derive the path integral which is of 
the form ~ e~ s , we need to "exponentiate" the bracket 
term in Eq. (|4.7[) and keep all the terms that are of or- 
der of At in the exponential. However, we cannot simply 
exponentiate the third term in the bracket because this 
term contains A which is of order of vAt. 

This is noted in Ref. [13], where the author derives the 
path integral for the Stratonovich convention, and cir- 
cumvents this difficulty by keeping the second order term 
in Ace in the exponential and replacing this term with its 
average value. Although the final expression is correct, 
that derivation might be inconsistent with the concept 
of path integral since that derivation is valid only in the 
mean-squared sense instead of for all paths, as required 
by the path integral. Here, we provide an alternative 
derivation that is valid for each path. First, we note that 
the last term in the bracket can be written as 
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dki 
2ngi 



a— (Axi - At fa 
9i 



+iii(Aa; i -At/ i )-ifef At _ 



"ft 



dki 

dki 
2ngi 



dk t 



+i|i(Ax i -Ai/ i )-ife?At 



(4.8) 



where the last line explicitly of order of At and can, therefore, be exponentiated without incurring any error to the 
first order in At. Returning to the conditional probability, we have 



P(Xiti\Xi-lti-l) 



+j|i(Ax i -At/ i )-ifc?At 



1 — a At f'i + laAtkig'i ) 



dki 
2ng l 

dkj e +«|i(Axi-At/ i +aAtg i 9<)-ifc?At-aAi/,; 



2ng l 
1 



At [ Agj 

5^4 "AT 



v / 2nAt gt 



(4.9) 
(4.10) 
(4.11) 



where the last line is valid to the first order in At. It should be noted that the Fokker-Planck equation, Eq. (|2.17l) . 
can also be derived using Eq. (|4. 1 1[) and the identity of Eq. (|4.3|) . This is done in the Sec. IIV Al Returning to Eq. 
P~4"|) , we have 



P(xftf\x t ) = 



dxi 



\j2itAtgx 



dx 



N-l 



1 



At 
^7 



\f2irAtg N ^ 1 \j2itAtg 



N 



Vx e 



-f*0 dt{^-±-^[d t x-f(x)+ a g(x)g>(x)Y+ a f>(x)) = f Xf 



Vx e 



(4.12) 
(4.13) 



with Vx = YliLi 



dxi 



v / 27rAtg i 



and the action given by 



S = 



dt 



it, 



2g{xf 



[d t x- f(x) + ag(x)g'(x)\ +af'(x) 



(4.14) 



where we have taken the formal limit of by letting N — > oo and At — > 0. Note the extra terms in the S coming from 
the Jacobian \dh/dxi\; they are needed in order to ensure that J dxfP(xftf\xo to) = 1- This can be demonstrated by 
explicit, but tedious, calculation for general a (see Appendix [B|. From Eq. (|4.14[) , it is clear that the Ito convention 
with a = is the simplest to deal with. Another particular useful form of the path integral is obtained using the 
Hubbard-Stratonovich transformation, which linearizes the quadratic term 



P(x f tf\x 1 ) 



Vy 



Vx exp - 



dt 



y 2 (t) - iy{t) [d t x - f(x) +ag(x) g'{x)] + af(x) }, (4.15) 



where the measure now is J Vy J** Vx — J . . . J ^ J c?xjv_i . . . J dx\. This result could, of course, also have 
been obtained directly by substituting ki = gi yi in Eq. (I4.10|l and taking the continuum limit. Note that in the 
discretized version of Eq. (|4.15p . y n is associated with time t* = at n + (1 + a) t n -i. This form of the path integral 
is closely related to the MSR formalism [32] to calculate response and correlation functions. This will be explored in 
Sec. ELS , 



It is interesting to see how the additional drift tern 
fi(x) in the Langevin Equation [Eq. (|2. 1 arises from 
the constraints that equilibrium statistical mechanics im- 
pose on the path integral formulation [35] . Thermal sys- 
tems must obey detailed balance which states that 

P(x f tf\x t )P eq (xq) = P{x t f \x f t )P eq (x f ). (4.16) 

The equilibrium distribution has the form P eq (x) = 
exp — l/37i(x)], and P{x$ tf\xf to) is the conditional prob- 
ability for the reversed path, i.e., for x{t) = x(—t). It 



turns out that the Stratonovich convention is the sim- 
plest for the discussion of time-reversal properties not 
only because it obeys the ordinary rule of differential 
calculus, but also because it has the property that the 
forward and backward paths are evaluated at the same 
points. We will employ the Stratonovich convention be- 
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low. First, we note that 

Peq(Xf) 



Peq(x ) 



exp-\pH(x f ) -pH(x )] 

d/3H{x) 



exp - 



dt (d t x) 



dx 



(4.17) 



and that P(xq tf\xf to) can be obtained simply by not- 
ing that the path associated with this distribution is the 
time-reversal path of P(xf tf\xg t ), which can be written 
as 



P(x tf\xft ) = T>y 



f-Tf 



Vx 



(4.18) 



x e 



■ fto dt l4 S 2 +nr[%z+/-i<7 9 ']+i/'] 



Now, using Eqs. l(i3B]) . (f4^T7|) . and P~T8"]) and comparing 
this term by term with exponential in Eq. (|4. 15[) [in the 
Stratonovich convention a = 1/2], we see that 



*»(*) 



8(3U(x) 
dx 



(4.19) 



f(x) = 2 9(x)g (x) ~ q£ ■ (4.20) 



The first term in Eq. (|4~20"1) is identical to Eq. ([2~2T]) in 
the Stratonovich interpretation. The second term is the 
standard frictional term, from which we identify the dis- 
sipation coefficient as T(x) = (3g 2 {x)/2, which is the Ein- 
stein relation. This derivation again demonstrates that 
equilibrium distribution is the only physics needed to fix 
f(x) for a given stochastic interpretation. 



A. Derivation of the Fokker-Planck Equation from 
the Path-integral 

In this subsection, we derive the Fokker-Planck equa- 
tion directly from the conditional probability Eq. (14. lip , 
thereby establishing the equivalence of the path integral 
formulation and the Fokker-Planck equation for general 
a. Let us rewrite the conditional probability, Eq. (|4.10[) . 
where we set x — Xi, t — ti, x' — Xi-i, t = and 
h = g l y l : 

P[xt\x't') 

-At [4 Wj+a/«'] 



2tt 
2tt 



* e ~AtAi+iyiAxi 



(4.21) 



where 



Ai[xi,Xi-i;yi] = y y\ + iyi [f% - ag^ + a/-, (4.22) 
Axi = x — x', and At = t — if . Our aim is to calculate 
P(x,t) = I dx'P(xt \x't')P(x',t'), (4.23) 



to first order in At. Expanding P(x' , t) 

P{x',t') = P(x-Axi,t') 

8 (At ) 2 8 2 

= P(x,t') A Xi g P (x,t') + ^L^ P{x ,t'), 



and putting this back to Eq. (14.23|) . we find that it can 
cast in the form 

P(x,t) - P(x.t') , . . , s n , . 8 , , N 

At = a(x)P(x,t')+f3(x)—P(x,t') 

+ j(x) ^P(x,t'), (4.24) 

where a(x) = limAt^o [1o{x) — 1] / At, f3{x) = 
-limAt^o 1\{x)/At, and j(x) = liniA^o 2^ (x)/(2 A*) . 
The main task is to evaluate integral of the form 

l m (x) = J ^ J dA Xi (Ax^Qix^Ax^y^e^, 

(4.25) 

where 



Q{x, Ax,; y{) 



-AtAi[x,x';yi] 



yi^L J^ Q{X Ax-y- 



(4.26) 

Axi=0 



where in the last line, we have Taylor expanded the func- 
tion Q(x, Ax^, yi). It is easy to see that 



8 m 

dA Xl {A Xl ) m e 1 ^' = 2ir(-i) m — S( yi ), (4.27) 



and therefore 



k-\-7n Qm-\-k 



, (4.28) 



y,=o 



^\ax>=o- Us " 



where Q^ k \x,0;yi) = 8 k Q(x, Ax l ;y)/8Ax k ^ 

ing Eq. (|4.28j) , it is straightforward to compute I m (x) to 
the first order in At; we obtain 

l (x) = 1 -f'(x) At 

+ (I -a) {[g'(x)} 2 + . g(x)g" (x)} At, (4.29) 

l^x) = [f(x)-(2-a)g(x)g'{x)}At, (4.30) 

l 2 (x) = g 2 (x)At, (4.31) 

with vanishing higher order terms, i.e. T n (x) = for 
n > 3. Therefore, we have 

a(x) = -f(x) + (l-a) [g'{xf + g(x)g" (x)](£.32) 
P{x) = -f(x) + (2-a)g(x)g'(x), (4.33) 

7 (x) = \g(x) 2 . (4.34) 

This is equivalent to the Mori expansion |3j. It is clear 
that with these coefficients Eq. (|4.24l) becomes Eq. (|2.17l) , 
the Fokker-Planck equation. 
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B. Correlation, Response functions, and 
Fluctuation-Dissipation Theorem 

One of the advantages of the particular form of the 
path integral in Eq. (|4.15p is that correlation and re- 
sponse functions can be computed conveniently from it. 
The average of any functional 0[x(t),y(t)] of x(t) and 
y(t) at fixed xq is given by 



(O[x{t),y{t)]) xo = Vy / VxO[x(t),y(t)]e-%4.35) 

J Jx 

In particular, the two-point correlation function is 



(*(*i)a:(t2)) ao = / Vy j Vxx{t t )x{t 2 ) e" 5 , (4.36) 

and the propagator function 

G{t 2 ,h) = {x{t 2 )[-iy{ti)]) X0 

= fvy f Vxxih^-tyit^e- 3 ^) 

J Jxq 

Physically, the propagator describes the response of the 
system to a delta perturbation. One of the nice features 
of the propagator function, which is useful in perturba- 
tive expansions, is that causality is automatically built- 
in, i.e. G(t, if) = if t < t'. To see this, we go back to the 
discretized form of the path integral and write G(t, t') as 



G{trn i in 



dx 



N ■ 



dxi 



dyN f dm 
2ir "J 2ir 

x x m [-i Vn ]e-^ AtA ^^ Ax \ (4.38) 



where Ai[xi, j/i] is defined in Eq. (|4.22[) . First, let 
us consider t n > t m ; each pair of the integrals J J dxi 
in Eq. (|4.38|) gives 1 for i > n. When integrating over 
x n , we make use the following identity 



J dx n Q\x n — \ , /\x ni yn]& 



+iy n Ax ri 



which gives zero when integrating y n . Thus, we have 
shown (x m [— iy n ])x = for all n > m. Now, sup- 
pose m — n, one can show that using the above iden- 
tity, (x„[-iy n ]} Xo = 1. Clearly, {x m [-iy n ]) Xo ^ 0, if 
to > n. Thus, we have shown how the path integral en- 
forces causality, i.e. G(t, if) = if t < t', and G(t, t') ^ 
if t > t' . However, there is a subtle point about the value 
of G(t, t) in the continuum limit, which has to be con- 
sistent with the a-convcntion. The simplest way do this 
is to note that since y n is really associated with time at 
i* = at n + (1 — a) t n —i, we have to evaluate 

G(t,t) = (x(t* n )[-iy n ]) X0 

= ( [ax n + (1 - a) x n -i](-Wn))x = a- 



Now, we specialize to a system near equilibrium, and 
we investigate how the path integral describes properties 
such as the Fluctuation-Dissipation Theorem. The equi- 
librium average of any function 0[x(t),y(t)] of x(t) and 
y(t) is defined as 



(0[x(t),y(t)}) 



eg 



Vy / VxO[x{t),y(t)}e- s P eq (x ). 



(4.40) 

Note that equilibrium averages are independent of a, pro- 
vided that we add the additional drift f\{x). When the 
system is under a time-dependent physical force h(t) 7 the 
total Hamiltonian is Ht = T~io(x) — x(t)h(t), so that 



f(x,t) = (1 - a)gg' - T(x) 
fo(x) = (1 - a)gg' - T(x) 
Therefore, we have 



dx 
dx 



= f Q (x)+T(x)h(t), 



6(x(t)) 



5h(t>) 



1 



h(t)=0 



2k B T 



{x(t){- ly {t')g[x{t')f}) 



eq 



k B T 

Xxx{t} t ) 



x(t)g[x(t')]g'[x(t')]) eq (4.41) 



We observe that the response Xxx(t,t') to a physical 
forces and the propagator G(t, t') defined in Eq. (|4.37|) are 
different, although they are proportional to each other 
for the case of uniform diffusion constant. In particular, 
there is an additional term arising from the normalization 
factor /' in the action and it is absent when the diffusion 
constant is spatially uniform. By integration by parts, 
the first term in the bracket can be evaluated to be 

{x{t){^y{t')g[x{t')f}) eq 

= ( x(t) {d t ,x(t') - f [x(t')] + ag[x{t')]g'[x{t')]}) eq . 

Therefore, the physical response function is 
1 



2k B T 



[x(t){d t ,x(t')-B{x(t')]}) eq , (4.42) 



where B(x) = gg' — T(x)dH.o/dx. Note that the physical 
response function is independent of a, as it should be; 
note also the a drift proportional to gg 1 arises from spa- 
tial varying diffusion constant. To proceed further, we 
note that as a consequence of the detailed balance con- 
dition, Eq. (|4.16[) . the equilibrium correlation function is 
symmetric with respect to exchange of t <-> if: 



(0 1 [x(t)]0 2 [x(t')}) eq = {0 1 [x(t r )}0 2 [x(t)]) 



eg 



Applying this to Eq. (I4.42|) and subtracting the results, 
we have 

[d t - d v ] (x(t)x(t')) eq = 2k B T [ X xx{t', t) - X xx(t, t')] . 
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Since the correlation function is time translational invari- 
ant, we must have d t > (x(t)x(t')) eq — —d t (x(t)x(t')) 
Thus, 

d t (x(t)x(t')) eq = ~k B T [Xxx(*,0 - Xxx(t',t)} . 

This is the Fluctuation-Dissipation Theorem. To put it 
in a more traditional form, we note that Xwx(t',t) = 
when t > t , and we can write 

xUt - = -^fdt (x(t)x(t')) eq 6(t - t>), 

where 9(t) is the Heaviside unit step function. The 
Fourier transform of the response function is given by 

(x(u)x(-u;)) eq = lmXxx{uj)i (4.43) 

which is of the form that is commonly quoted in the 
literature. 



C. Perturbation Theory 

One of the advantages of the path integral formulation 
of stochastic dynamics is that it is by construction a field 
theory that facilitates systematic perturbative calcula- 
tion of correlation functions. In particular, for systems 
with state-dependent dissipative coefficients, the result- 
ing Langevin equation is generally nonlinear, and per- 
turbation theory is a convenient way to derive the mode- 
coupling theory [36| . Thus, in this subsection, we set 
up the perturbation theory for a systematic calculation 
of correlation and response functions in the deviation of 
the diffusion coefficient from spatial uniformity. First, 
we need to set up the generating functional. Note that 
Peq(xo) satisfies 



P eq {x)= J dx s P(xt\x s t s )P eq (x s ), (4.44) 

which implies that 

P eq (x)= lim P(xt\x s t s ). (4.45) 

t s —*- — oo 

Thus, the equilibrium averages can be written as 



(0[x(t),y(t)]) eq = I Vy I VxO[x(t),y(t)]e- s P eq (x Q ) 



Vy I VxO[x(t),y(t)]e- s , 



where in the last line, the limit of the time integration 
in the action S is extended to — oo to oo. This allows us 
to define the generating function for equilibrium averages 
by 



Z[F,F] = / Vy / Vx e~ s+i dt[x{t)F{t) ^ lv{t) ^ {t)] . 




= +r 7fc x 3 x 




2-2 



FIG. 3: The two vertices corresponding to Eq. (|4.49|l for a 
particle diffusing in a spatially varying diffusion coefficient 



given by D{x) 
potential. 



Dq(1 + yx ), and confined in a harmonic 



The correlation functions and the propagator are simply 
functional derivatives of Z. This sets up the MSR per- 
turbation scheme [32j that allows the immediate applica- 
tion of all of the powerful techniques of field theory, in- 
cluding the renormalization group, to nonlinear stochas- 
tic problems. It is customary to introduce the variables 
x(t) = —iy(t). Note that in the perturbation expansion, 
all the a dependent terms cancel provided that we use 
(x(t)x(t))o = a (see Appendix |C|) . Therefore, it is con- 
venient to use a = at the outset. 

As an informative model calculation, we explore the 
problem in which a particle diffusing with D(x) — 
£*o(l+7^ 2 ) confined in a harmonic potential, TL = 
kx 2 /2. If the confining potential is turned off, this prob- 
lem is exactly solvable, as shown in Sec. lIII Cl It can also 
be solved exactly when 7 = but not when 7^0, and 
a perturbative expansion in 7 is useful. The goal of this 
exercise is to compute the propagator (xx) and the cor- 
relation function (xx) separately and to check that the 
Fluctuation-Dissipation theorem is satisfied. For simplic- 
ity, we work in the Ito convention (see Appendix [C] for 
general a), and set ks = 1- According to the formalism, 
we have 



f{x) = gix)g\x) 

= -T ik-2yT)x 



TV \ m 



r kyx 3 



(4.47) 



It should be pointed out that from Eq. (|4.47p one might 
at first sight conclude that there is a broken-symmetry 
state, when k < 2yT, with (x) =/= 0. But we know that 
this cannot happen because the stationary distribution is 
in fact the Boltzmann distribution. Therefore, one could 
get the wrong physics if one only looks at "classical" tra- 
jectory, i.e. solution to dt,x = fix), which maximizes 
the action S in the Ito convention. This shows again the 
importance of noise in these problems. 

The unperturbed and perturbing actions are 

So = J dt [-D a x 2 + x id t x + T k'x)] (4.48) 

5 , = r 07 /*[*A-r,V], (4.4 9 ) 



(4.46) 



where k' 



27T. Introducing the state vector x c 
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(b) 




(c) 



Id) 



FIG. 4: Diagrams that contribute to the self-energy S(oj). 
Note that diagram d is identically zero. 



(x, x), we can write Sq as 
1 f dio 



2 / 2tt 



where 



3* M 



-2D -iuj + T Q k' 
iuj + T k' 



(4.50) 



and thus 



^aM^J - 1 2_D , 

y -iw+r fe' piw+rw? / 

from which we can read off the bare propagator and the 
zeroth-order correlation function: 



(x (v)x (o/)) 
{x(lo)x{lo'))q 







1 



—iuj + Tok 

\-iuj + r k'\ 2 



The interacting Si consists of two vertices that are de- 
picted in Fig. [3J To second order in 7 , the inverse of the 
propagator, G~ 1 {uj), can be written in frequency space 
as 



G-\lu) = — 1 u> + T (fc + 7 T) + E(w). 



(4.51) 



The self-energy are computed from the diagrams 

listed in Fig. 5] and it is given by 



= 6fc(r o7 ) 2 [2TA(u) - 3kB(u>)\ 



where 



A{uj) = J — ^ J G {uj - LUl ~ UJ 2 )G Q {uJl)C a {uJ2), 

du>i f du>2 



B{u>) 



2vr 



2tt 



G (uj -u>i- ll> 2 )C {uji)C (u>2), 



where Cq{lo) = {x{lo)x{—uj))q is the zero-order correla- 
tion function. After some algebra, we obtain 

= -iuj + r fc [ 1 + jT/k - 6 ( 7 T/fc) V^/r fc) 1 , 



(b) 




(c) 



(d) 



FIG. 5: Diagrams that contribute to the noise D(uu) 



where /i(s) = (—is + 3) _1 . 

Next, we turn to the correlation function C(t, t') = 
(x(t)x(t ) ), which can be written in the form: C(u/) = 
2D(w) |G(w)| , with D(u>) are computed from diagrams 
listed in Fig. [3] and it is given by 

D(u) = D Q fl + ^-j+kT 2 ol 2 [3kE(uj)-8TH(uj)} 



E(u>) 
H(u>) 



du>i 
"2V 
dw\ 
"2V 



du 2 

~27T~ 

~2^ 



C (uj -coi- u 2 )C (ui)C (u; 2 ), 



X [Gq(ijJ — W\ — UJ2) + Gq(uJi + U)2 — U}) ] . 

After some algebra, we find 

D{u) = 2D [ 1 + ( 7 T/fc) - 6( 7 T/fc) 2 Re/x(w/r fc) ] , 
and we evaluate the correlation function, 



T 
I 

T 
I 



3Q!2 — Oil 



a 2 1 



k 



3 a\ 



,(4.52) 



Q 2 



with decay rates 

2 

r fc 

r fc 
2 



4 + 7 T/fc - v/4-4 7 T/fc + 25 ( 7 T/fc) 2 
1 + 7 T/fc - 3 ( 7 T/fc) 2 j , for 7 T/fc < 1, 



4 + 7 T/fc + ^ 4 - 4 7 T/fc + 25 ( 7 T/fc) 2 
3r /c [l + ( 7 T/fc) 2 ] , for 7 T/fc < 1. 



Note that there are now two decaying modes with a fast 
mode 012 and a slow mode ol\ in the system in contrast 
to the case with with uniform diffusion. Note also that 
(x 2 (0)) = T/k as it should be. If we did not put in the 
extra drift term gg' in f(x), this relation would not hold. 
In fact, it would have been (x 2 (0)) = T/(k - 7 T), which 
violates the equipartition theorem. This is yet another 
demonstration that this extra drift term gg' is needed to 
ensure the correct thermodynamic properties. In Fig. [5J 
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FIG. 6: A plot of the correlation function of x, C(i,0), as 
a function of t for yT/k = 0.15. The solid line represents 
C(t,0) as given by Eq. (|4.52[1 . which is calculated from the 
perturbation theory to second order in jT/k. The data points 
are obtained from numerical simulation of the corresponding 
Langevin equation. Clearly, the result from the perturbation 
theory agree excellently with the simulation. The dashed 
line represents the bare (zeroth order) correlation function, 
which has a decaying rate that is different from the case if 
the spatial- varying dissipative coefficient. 



we plot the correlation function in Eq. (|4.52p and the nu- 
merical simulation of the Langevin equation describing 
this system for ^T/k — 0.15. Clearly, the second order 
perturbation theory agrees very well with the simulation. 
Note, however, that when jT/k ~ 1, a\ becomes nega- 
tive, signalling the breakdown of perturbation theory. 

Finally, we demonstrate FDT to second order in per- 
turbation theory. The physical response function is given 
by 

Xz*(tJ) = T (x(t)x(t')[l+ 7 x 2 (t')]e- Sl ) Q 
= r a G(t,t')+T ol G o (t,t')C o (0) 
- r o7 ( a ;(t)x(t'^ 2 (i')5/) , 

which corresponds to the diagrams in Fig. [7] We find 

XM = T Q G(uj) [1 + yT/k - 6 ( 7 T/k) 2 fi(cj/T k)] . 

(4.53) 

Note that this clearly shows that the physical response 
function and the propagator are different. Taking the 
imaginary part of Eq. (|4. 53[) , it can be easily verified that 
the Fluctuation-dissipative theorem Eq. (|4.43)1 is indeed 
satisfied to second order in perturbation theory. 



V. N-COMPONENTS LANGEVIN EQUATION 

Many physical problems involve more than one vari- 
able and some of the issues we have addressed so far 
may not apply to higher-dimensional systems. For ex- 
ample, the drift term f{x) in 1-D can always be written 
as a derivative of another function, i.e. 1-D systems are 
conservative, however, for higher dimensional systems, 
this may not be true. A complete analysis of higher di- 
mensional systems requires a separate publication. Here, 



(a) 



(b) 



) ^ (d) 



(e) 



FIG. 7: Diagrams that contribute to the physical response 
function 



we briefly discuss the Fokker-Planck equation and the 
path integral in the a-convention for a multidimensional 
Langevin equation of the form 

d t Xi = fi(xi, ■ ■ ■ ,x N ) +gij(xi,- ■ ■ ,x N )rjj(t), (5.1) 

where f]i(t) are the noises, with zero mean and correlation 
given by 



(r h (t)r )j (t')) = 6 ij 6(t-t'). 



(5.2) 



In Eq. (|5.1|) and the following, Einstein summation is 
assumed. We focus on the case where the system is 
near thermal equilibrium and address, as we did in the 
ID case, how the Boltzmann distribution determines the 
form of fi and gij in the a-convention. The Fokker- 
Planck equation corresponds to Eq. (|5.1[) can be derived 
following the same procedure as outlined in Sec. Ill CI In 
the a-convention, we find 



dV 
~dt 



d 

dxi 



f , 9 9U 
J i + " "5 — 9ki 



(5.3) 

where 7 3 [{x},t] is the joint probability distribution of x$ 
at time t. If there are only dissipative terms and no 
reactive terms in the constraint that V reach a long- 
time state of thermal equilibrium value proportional to 
exp — j3H requires that 



1 d 



dgij 

a^-gij 



dxi 



in order that Eq. (|5 . 3|) reduce to 



dV 
dt 



_d_ 



2#ii 9ji 



dxj dxj 



i on 

2^ 9ik9 ' k dx~' 
(5.4) 

(5.5) 



dV 



with the steady-state solution V cx exp— PH. 
The diffusion matrix is defined as 

D ij({x}) = -^gugn = k B TT i:j ({x}), 



(5.6) 
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where Tij({x}) is the matrix of dissipative coefficients. 
Thus 



fi({x}) 



dD t3 ({x}) 
dxj 



dgij 



3xa ' 



(5.7) 



be symmetric in which case, it is simply the square root 
of /c B 7T l:) 

To derive the path integral, we first discretize Eq. ([57 

as 



Note that since the diffusion matrix Dij({x}) is a sym- 
metric with respect to i «-> j, it has only N(N + l)/2 inde- 
pendent entries, and we may impose iV 2 — N(N + l)/2 = 
N(N — l)/2 constraints on gij({x}) without sacrificing 
the physical content. We could, for example chose to 



ct ] = xt^ + Mf^ + gP I *»&(*)> (5 



and introduce 



{4 n) },{4 n - 1) }]-[^ 1 



where f- n> = fi({ax\ n ' + (1 — a)x\ n 1> }) and g^' — gij({ax\ n> + (1 — a)x\ n 1> }). Following the basic steps as outlined 
in Sec. HVi we can write the conditional probability as 



„(«-i)l 



» 



At/ 



(n) 



tn 



dsr}i(s), 



(5.9) 



tn-l 



P({4 n) }in|{4 n_1) }*n-0 = det 



dx 



(n) 



I N 

n*[^({4 w) >.{*ir" l) }) 



\i=l 



(5.10) 



Taking the derivative of hi explicitly, we find 



dh 



dx 



where we have defined the matrix Mik by 



\a - a + 5 ^ , d 9lm -i 

Mlk = aM dx- k +a ^T 9 - 



x (n) _ fn-1) _ Mf (n) 



V k OXk 

Using the identity dctX = cxpTrX, the determinant can be evaluated to give 



det 



dhi 



dx 



(n) 



1 



det gij 



1-M U + - (MuMkk - MikMu) 



(5.11) 



(5.12) 



Therefore, the conditional probability can be written as 

N 



P({4 n) }*n|{4" _1) }«n-l) = n 



dki 



det g^ 



l-M u + - [MuMkk - MikMki) 



(5-13) 

where we have only kept terms up to order At. Following the similar procedure leading to Eq. (|4.11|) for the 1-D case, 
we find 



N 

^({4 n3 }*ni{4 n " i3 }*n-o = n/^ 



dkj 



det 



9a 



1 - aAtdifi + aik m digi m At + — At {d k gi m dig km - digi m d k gkm) 



x e 



%k i gT^(Ax j -At} i )-\k^ At 



(5.14) 



f ^Mt. p+ l yk{Ax k -f k At+ag k jdigi j At)-±gi k g kj y l y j -aAtdifi- 2 §-(d k gi,- n dig km -digi m d k g km )At 



2tt 

where in the last line, we have exponentiate terms in the bracket and substituted yi — gijkj. In the continuum limit, 
we have 



P({4 /) }*/K4° ) }*o)= L T>x k Vy k e- S , 

Jx[ J 

1 a 2 
S = I dt - gikgjkyiU] - Wk ( d t x k - fk + a gkjdigij ) + a dji + — ( d k gi m dig km - digi m d k gkm) 



(5.15) 
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Note that there is an extra term proportional to a 2 . This term is identically zero for 1-D system. If we start with 
Eq. (|5.14[) and follow the procedure as outlined in Sec. IIV Al we can show that the path integral is equivalent to the 
Fokker-Planck equation in Eq. (|5.3j) . 



VI. CONCLUSION 



In this paper, we have examined a thermodynamically 
consistent Langevin formulation of the Brownian mo- 
tion with a diffusion coefficient that depends on space. 
We argue, in particular, that the requirement that the 
Boltzmann distribution be reached in equilibrium deter- 
mines the interpretation of stochastic integrals arising 
from multiplicative noise in the Langevin equation. We 
hope that this paper clarifies some of the confusion over 
these stochastic issues that have persisted for some time. 
We have also constructed path integral representations 
of the Langevin equations with multiplicative noise, and 
we used this representation as a starting point for the 
development of a systematic perturbation theory. Such a 
formulation can be employed to treat nonlinear stochastic 
equations arising from a variety of problems. Future work 
includes generalizing this formulism to "fields" and ex- 
amines how state-dependent dissipative coefficients may 
give rise to long-time tails and corrections to scaling in 
dynamic critical phenomena. Of course, one of the most 
interesting open questions is whether there is an equiva- 
lent criteria for systems that are driven far from equilib- 
rium. 
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APPENDIX A: CONNECTION BETWEEN a 

AND 6»(0) 



In this Appendix, we outline the connection between 
a and the Heaviside unit step function, 9(t) evaluated at 
t = 0. For simplicity, we set f(x) = in Eq. (|2~T0|) : 



Using the a-convention rule Eq. 
Ax(t) = x(t + At) - x(t) 



we have for 



/t+At 
dsrj(s) 



/t+At 
dsr)(s) 

/t+At 
ds-q(s) H 

/t+At 
dsr)(s) 

/t+At pt+At 
ds / ds'ri(s)ri(s l ) + ■ ■ ■ 

Therefore, the average Ax(t) over the noise is 

pt+At pt+At 

(Ax(t)) = agg' ds ds'5(s-s'). (A2) 



The integral 



ds' S(s — s') 



s > t 

1 s < t 



= 0(t-s), 



(A3) 
(A4) 



defines the Heaviside unit step function, 9(t), and Eq. 
121) becomes 

pt+At 

(Ax(t)) = agg' / ds [6(t + At - s) - 6{t - s)] 



(A5) 



agg' At. 



On the other hand, we can directly integrate Eq. (|A1[) to 
obtain 

/t+At 
dsg[x(s)]r)(s). (A6) 

Expanding <7[x(s)] as 

g[x(s)} = g[x(t) + Ax(s)} = g[x(t)\ + g'[x(t)]Ax(s) - 
we find 



(A7) 



pt+At pt+At ps 

Ax(t) = g v( s )+gg' ds / ds'r}(s)r](s'). 

Jt Jt Jt 

. (A8) 

Note that the upper limit of integration for s' is different 
from that in Eq. (|A2[) . Therefore, we find 



(Ax(t))=6(0)gg'At. 



(A9) 



d t x = g(x)r)(t). 



(Al) 



Comparing this with Eq. (|A5|) . we conclude that 6(0) 
a. 
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APPENDIX B: EQUIVALEN CE B ETWEEN EQ. 
ff47fT> AND EQ. (|47TTj) 

In this Appendix, we demonstrate that the conditional 
probability given in Eq. (|4.7|) and its exponentiated form, 
Eq. (|4.1ip . are indeed equivalent. Note that the latter ex- 
pression involves a subtle step which is required for the 
construction of the path integral in Sec. HVl Therefore, 
it is crucial to confirm Eq. (|4.11j) is correct at least to 
order of At. We have already checked that the correct 



Fokker-Planck Equation, Eq. (|2.17p . can be derived from 
Eq. (|4.11[) in Sec. IIV At Here, we check that the normal- 
ization condition: 

J dxiP(xiti\xi-xti-i) = 1, (Bl) 

for general a is satisfied by Eq. (|4.11[) . First, let us check 
Eq. (|B1[) is true for Eq. (|4.7[) . For simplicity, we set 
f(x) = 0. We have 



dx i P{xit i \x i _xti_ l ) 



dA Xi e -%myr 



\J1-KAtgi 



l-a^Ax, 
9% 





" dw 


-J 





1 - a^-VAtw 
9i 



(B2) 



where we have made a change of variable: Axi — \ Atw. Remembering that gi = g[xi-% + a\ Atw] and expanding 
them in Eq. (|B2|) in power of At, we have 



dxiP{xit i \x i -xt i -x) 



dw -™1 
e ^ 



2tt 5 



1 - 



2ag' 



( 3a 2 g' 2 3a 2 g " 

[~ 2 27 



Atw 



At', 



ag 



Atw J 



a 2 g" 7a 2 g 

~W 3 V 



a 2 g' 2 



2g 6 



Atw b 



= l + 0[At 



3/2] 



(B3) 
(B4) 



where g = g(xi~i). Thus, Eq. (|4. 7[) indeed satisfies the normalization condition, which is hardly surprising since it 
must be true by construction. Now, let us check the exponentiated from, Eq. (|4. 1 1[) . We have 



dx^P {xiti I Xi—\t 

dw 



dAa 



y/2TrAtgi 



2ng 



■e 2 f 



i-l - 

J V^TT 

-, 1 2 a 2a g' HT7 
1 a n V At w 

2 y 



(la 2 g 12 ia 2 g" 



dw -fA\vm +a ^) 



g 



\ 2. 9 2 



Atw 



2g 6 



l + 0[At 



3/2l 



2.9 



Atw^ 



ag 



Atw 6 



a 2 g" 7a 2 g 



1„I2 



2<? 3 2 3 4 



(B5) 
(B6) 



Thus, Eq. (|4.1ip also satisfies the normalization condi- 
tion. We note in passing that although the expansions, 
Eqs. (|B3|) and (|B5[) are different, they both give one at 
the end result to the lowest order and the next order term 
is of the order At z / 2 . 



APPENDIX C: ALPHA-DEPENDENT 
PERTURBATION THEORY FOR THE MODEL 
SYSTEM WITH D(x) = D {1 + -yx 2 ) 



In this Appendix, we carry out a first order perturba- 
tion calculation for general a of the model system studied 
in Sec. IIV Cl in which D(x) = Dq(1+^x 2 ) in order to clar- 



ify some subtle issues associated with the a-convention. 
It is straightforward to work out the action in Eq. (|4. 1 5[) . 
Up to an irrelevant constant, we have 



dt [-D Q x 2 + x (d t x + T kx)] 



50 = 

51 = T 7 J dt 

— 3 a — x 
T 



2(2a - l)xx + — x 3 x 



2-2 
X X 



Note that there are more diagrams to evaluate than 
there are for a = 0. Consider first the propagator 
G(t,t') = (x(t)x(t')) , which can be written as G~ 1 (w) = 
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(a) 

www- 



(b) 




(c) 




FIG. 8: Diagrams that contribute to the self-energy E(o»). 



Gq 1 (uj) + S(o>), where the diagrams for the self-energy 
are displayed in Fig. [8] Note that the closed loop dia- 
gram c in Fig. [5] contains Go(t = 0) which must be set to 
a as explained in Sec. IIVBI We find 

E(w) = 2{2a-l)D o -f + 3D Ql -4D ol G Q {0) 
= D -f, 

which agrees with the calculation for a = 0. Note that 
the final result is independent of a as it should be. To 
first order, the noise D(u>) renormalizes exactly the same 



way as in the a = calculation. However, the physical 
response function is different. It is given by Eq. (|4.41[) 
with an extra a dependent term: 



Xxx{t,t') = T Q {x(t)x(t')[l + 1 x 2 (t')]e- s ') Q 



a 
T 



(x(t)[2 7 D x(tf)]e- s % 



7 T 



To \1 + ^ G(t, t') + 2aT ol C Q (t, t') 
2ar O 7G,(M')+0(7 2 ) 
k 



To 1 



G(t,t'). 



Without the cancellation of the a dependent terms, 
Xxx{t,t') would not have been causal. 
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